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The equilibrium structures, harmonic vibrational frequencies, quadratic force fields, dipole moments, and IR 
intensities of several triatomic bromine compounds of known or potential importance in stratospheric ozone 
depletion chemistry have been determined using the CCSD(T) electron correlation method in conjunction 
with a basis set of triple zeta double polarized (TZ2P) quality. Specifically, the molecules included in the 
present study are HOBr, HBrO, FOBr, FBrO, BrNO, BrON, B^O, BrBrO, BrCN, BrNC, ClOBr, CIBrO, and 
BrClO. Very accurate isomeric energy differences have also been determined at the CCSD(T) level with 
atomic natural orbital basis sets that include through g-type functions. In most cases, the isomer with a 
normal neutral Lewis dot structure is the lowest energy form, with the single exception that FBrO is predicted 
to be 1 1.1 kcal/mol (0 K) lower in energy than FOBr. In all cases, however, the hypervalent isomer is more 
stable relative to the isomer with a normal Lewis dot structure as compared to the chlorine analogs. Consistent 
with this observation, the energy of the last three molecules given above increases in the order ClOBr < 
CIBrO < BrClO. The CCSD(T)/TZ2P geometries and vibrational frequencies are in good agreement with 
the available experimental data. Heats of formation are determined for all species using a combination of 
theoretical isomeric, homodesmic, and isodesmic reaction energies. The accuracy of these quantities is 
ultimately dependent on the reliability of the experimental heat of formation of HOBr. 


Introduction 

There is growing interest in developing a more complete 
understanding of stratospheric chemisty, especially halogen and 
nitrogen oxide chemistry that leads to ozone depletion. As such, 
there has been considerable work over approximately the last 
30 years aimed at better characterizing chlorine, fluorine, and 
nitrogen oxide molecules that are prevalent in the stratosphere. 
Over the last few years this laboratory has been involved in the 
ab initio characterization of several chlorine, fluorine, and 
nitrogen oxide species (e.g., see refs 1—5 and references therein). 
Although bromine is a minor constituent of the stratsophere, it 
is known to have a greater ozone depletion potential than 
chlorine (for example, see refs 6 and 7, and references therein). 
Until fairly recently, however, bromine compounds have at- 
tracted much less interest, and the present ab initio study is 
designed to provide a better characterization of triatomic 
bromine-containing oxide compounds that are of potential or 
known importance in stratospheric chemistry. Specifically, the 
HOBr, FOBr, BrNO, BrON, Br^O (i.e., BrOBr), BrCN, BrNC, 
and ClOBr molecules and the hypervalent HBrO, FBrO, BrBrO, 
CIBrO, and BrClO compounds have been studied via ab initio 
calculations in the present investigation. 

Nitrosyl bromide, BrNO, is probably the best experimentally 
characterized of the molecules included in this study. The 
vibrational spectrum, 8 - 9 microwave spectrum and molecular 
structure, 10-12 electronic spectrum, 33 and heat of formation 34 of 
BrNO have all been the subject of previous investigations. At 
this point, the gas-phase fundamental vibrational frequencies 
and a vibrationally averaged molecular structure are both well 
established and therefore BrNO will serve as a benchmark for 
evaluating the reliability of ab initio geometries and frequencies 
for bromine oxide molecules that are obtained with the singles 
and doubles coupled-cluster method that includes a perturba- 
tional estimate of the effects of connected triple excitations. 
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CCSD(T), 15 in conjunction with a one-particle basis set of triple 
zeta double polarized (TZ2P) quality. The heat of formation 
of BrNO is reported to be well established, 34 and this will enable 
us to determine accurate heats of formation for other bromine 
compounds using a combination of isomerization and isodesmic 
ab initio reaction energies together with some experimental data 
on chlorine compounds. We are aware of only one previous 
ab initio study on BrNO in which the effects of electron 
correlation were explicitly included. Meredith, Quelch, and 
Schaefer 16 studied both the BrNO and BrON isomers at the 
singles and doubles configuration interaction (CISD) level of 
theory using a double zeta polarized basis set. Their CISD/ 
DZP equilibrium geometry was in modest agreement with the 
experimental structure, with the theoretical Br-N bond distance 
too short by ^0.04 A. In the earlier studies of chlorine oxide 
species 1-5 the CCSD(T) level of theory has proven to yield 
equilibrium geometries and vibrational frequencies that are 
generally more accurate than those at the CISD level of theory . 17 

The HOBr and BrCN molecules are the next most studied 
bromine triatomics included in this study. These are the only 
other species for which vibrationally averaged structures 18-20 
and all fundamental vibrational frequencies 20 * 21 have been 
experimentally observed or deduced. The heat of formation of 
HOBr has also been examined in several recent investigations, 22-24 
but the variation among these values is rather large. This is 
discussed in more detail later. We are aware of only one 
previously published ab initio study of HOBr in which the 
effects of dynamic electron correlation were included. In this 
study, McGrath and Rowland 22 determine a A//£ 300 value of 
-14.2 kcal/mol using G2 theory, although they do not specify 
the one-particle basis set that was used for bromine. Further- 
more, G2 theory 25 was not designed for molecules containing 
third-row atoms (such as bromine) and therefore this value must 
be viewed with caution. For BrCN, we are aware of only one 
previous ab initio study in which the effects of electron 
correlation were explicitly included. Kello and Sadlej 26 exam- 
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ined the electrical properties of the cyanogen halides (F, Cl, 
Br, and I) using large basis sets, sophisticated coupled-cluster 
electron correlation methods, and experimental geometries. They 
showed that relativistic effects do not become significant until 
ICN. 

Some of the fundamental frequencies of B^O, BrBrO, BrNC, 
BrCIO, and ClBrO have been observed experimentally 2127 and 
these are discussed further below. Properties determined in the 
present investigation include equilibrium geometries, dipole 
moments, harmonic vibrational frequencies, infrared (IR) in- 
tensities, isomerization energies, and heats of formation. The 
ab initio methodology used in this study is discussed in the 
next section and the results and discussion are presented in the 
following section. Conclusions are presented in the final section. 

Computational Methods 

The equilibrium geometries were determined with a TZ2P 
basis set at the CCSD(T) level of theory. The TZ2P basis set 
used for the C, N, O, and F atoms consist of Dunning’s 28 5s3p 
contraction of Huzinaga’s 29 10s6p primitive sets, and the H 3s/ 
5s set comes from the same source. The polarization orbital 
exponents (two p functions for H and two d functions for C, N, 
O, and F) are taken from Dunning. 30 The Cl TZ2P basis set 
consists of McLean and Chandler’s 31 6s5p contraction of 
Huzinaga’s 12s9p primitive set augmented with two sets of 
polarization functions with exponents, Od — 1.072 and 0.357. 
The Br TZ2P basis set is composed of a 6s5p2d contraction of 
a 17sl3p6d primitive set as given by Schaefer, Huber, and 
Ahlrichs. 32 The orbital exponents of the two d polarization 
functions are ctd — 0.674 and 0.225. All six components of 
the Cartesian d functions were included in the basis sets. 
Coupled-cluster analytical gradient methods 33 34 were used to 
locate equilibrium structures, while quadratic force constants, 
harmonic frequencies, and IR intensities were determined by 
finite differences of analytical gradients. The dipole moment 
was determined as the derivative of the energy with respect to 
an external electric field. In all coupled-cluster calculations, 
the C, N, O, and F ls-like core molecular orbitals, the Cl ls- 
and 2sp-like core molecular orbitals, and the Br Is-, 2sp-, and 
3spd-like core molecular orbitals were constrained to be doubly 
occupied in all configurations (i.e., the frozen core approxima- 
tion was used). In addition, the C, N, O, F, Cl, and Br Is core- 
counterpart virtual molecular orbital was deleted from the TZ2P 
basis coupled-cluster calculations. 

In order to determine accurate isomeric energy differences 
and heats of formation, reaction energies have been evaluated 
at the MP2, CCSD, and CCSD(T) levels of theory using very 
large atomic natural orbital (ANO) basis sets. 35 The ANO basis 
set for Br was recently developed by Bauschlicher, 36 while the 
ANO basis sets for the other atoms have been described in detail 
previously. 35 37 The basis set denoted ANOl consists of 6s5p3d, 
5s4p2d, 4s3p2d, and 4s2p ANOs on Br, Cl, (C, N, O, and F), 
and H, respectively, while the AN02 basis set is composed of 
6s5p3dlf, 5s4p2dlf, 4s3p2dlf, and 4s2pld ANOs on Br, Cl, 
(C, N, O, and F), and H. The AN03 basis set consists of 
7s6p4d2f, 6s5p3d2f, 5s4p3d2f, and 4s3p2d ANOs on Br, Cl, 
(C, N, O, and F), and H, respectively, while the AN04 basis 
set is composed of 7s6p4d2flg, 6s5p3d2flg, 5s4p3d2flg, and 
4s3p2dlf ANOs on Br, Cl, (C, N, O, and F), and H. For the 
ANO basis sets, only the spherical harmonic components of 
the d-, f-, and g-type functions were included. 

The coupled-cluster geometry optimizations were performed 
with the TITAN 38 program system. The MP2 and coupled- 
cluster single-point energies were performed with the TITAN 
coupled-cluster programs interfaced to the SEWARD 39 integral 


program, and the SWEDEN 40 self-consistent field and trans- 
formation programs. All calculations were performed on either 
the Computational Chemistry Branch’s IBM RS6000 590 
computers or the NASA Ames Central Computer Facility’s Cray 
C90. 

Results and Discussion 

A. Equilibrium Structures and Dipole Moments. The 

CCSD(T)/TZ2P equilibrium structures, rotational constants, 
dipole moments, and 5T\ diagnostics 41 are presented in 
Table 1. Experimental data for HOBr, BrNO. and BrCN are 
also given for comparison. The 57} diagnostic values for the 
species with normal neutral Lewis dot structures (i.e., HOBr, 
FOBr, BrNO, B^O, BrCN, and ClOBr) are all less than 0.020 
except for BrNO, which has a value of 0.023. The electronic 
structure of nitrosyl halides is somewhat more difficult to 
describe (e.g., see ref 4, and references therein), but a 5T\ value 
of 0.023 is still not very large. The two molecules possessing 
normal Lewis dot structures, but where this electron configu- 
ration has partial charges on two atoms (Br— N + — C' and Br— 
0+-N"), have larger 5T\ diagnostic values than their respective 
isomers possessing normal neutral Lewis dot structures, espe- 
cially BrON. It is noteworthy that for BrNC and BrON, the 
electron configuration having a normal charged Lewis dot 
structure does not seem to be the dominant configuration, at 
least based on Mulliken population analyses presented later. The 
5T\ diagnostic values of the hypervalent species (HBrO, FBrO, 
BrBrO, ClBrO, and BrCIO) are all significantly larger than the 
value obtained for the respective isomers possessing normal 
neutral Lewis dot structures. This is to be expected since it is 
known that the electronic structure of hypervalent species is 
more difficult to describe properly, but the important point is 
that the 5T\ diagnostic values of all the bromine species included 
in this study are well within the range where the CCSD(T) 
method is known to yield reliable geometries, vibrational 
frequencies, dipole moments, and other properties (e.g., see ref 
17, and references therein). 

As indicated previously, experimental geometries are known 
only for HOBr, BrNO, and BrCN, and these are all vibrationally 
averaged structures of one sort or another. Nonetheless, the 
agreement between the CCSD(T)/TZ2P equilibrium structures 
and the experimental vibrationally averaged structures is gener- 
ally good, with the ab initio bond distances being too long. The 
agreement for the HOBr and BrNO bond angles is excellent. 
In all cases, the Br— X bond distance exhibits the largest 
difference between theory and experiment, which is largest for 
the Br— N bond in BrNO (0.033 A). There is little experience 
on which to judge the error in the CCSD(T)/TZ2P values for 
Br— X bond distances, but the differences for the O— H, N— O, 
and C— N bond distances are all as expected based on the 
correlation treatment and one-particle basis set that was used 
(see ref 17 and references therein for examples; note that these 
errors will decrease as better one-particle basis sets are used 
with the CCSD(T) method). Our previous studies 1-5 on Cl— X 
bond distances, however, also showed similar behavior (i.e., 
that the difference between CCSD(T)/TZ2P and experiment is 
somewhat larger for Cl— X bonds than found for bonds between 
first row atoms) and therefore it is expected that the errors for 
the CCSD(T)/TZ2P Br-O, Br-N and Br-C bonds in Table 1 
are typical. In fact, the errors in CCSD(T)/TZ2P Cl— O, Cl— 
N, and Cl— C bonds are quite similar to those observed here 
for the analogous bromine bonds. Based on these comparisons, 
the theoretically predicted structures contained in Table 1 should 
greatly aid in the analysis of future experiments. It should be 
borne in mind, however, that also based on earlier studies of 
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TABLE 1: Total Energies (E h ), Equilibrium Structures, 

Rotational Constants (MHz), and Dipole Moments (D) of 

_ . . ai l 1 DoAminoS 


HOBr 


FOBr 


BrNO 


Br 2 0 


BrCN 


ClOBr 


BrCIO 


parameter 

CCSD(T) 

expt* 

molecule 

parameter 

CCSD(T) 

E 

0.220646 


HBrO 

E 

0.122600 

V? 

0.012 



7? 

0.023 


1.49 

1.38 


V 

4.00 

^HO 

0.964 

0.961 


r Hbt 

1.458 

/"OBr 

1.853 

1.834 


fBrO 

1.731 

ZHOBr 

102.3 

102.3 


ZHBrO 

106.1 

A c 

609990 

613677 


A e 

261720 

Be 

10364 

10577 


Be 

12526 

C e 

10191 

10383 


Ce 

11954 

E 

0.208952 


FBrO 

E 

0.216515 

7? 

0.018 



5? 

0.029 

(A 

1.58 



fi 

2.68 

rpO 

1.464 



r FBr 

1.841 


1.846 



'‘BrO 

1.661 

ZFOBr 

106.3 



ZFBrO 

109.0 

Ac 

44695 



Ac 

20859 

Be 

4036 



Be 

6951 

Ce 

3702 



Ce 

5214 

E 

0.249567 


BrON 

E 

0.213607 

*7? 

0.023 



v? 

0.036 

pL 

2.06 

1.80 


u 

2.78 

rBr N 

2.173 

2.140 


rbto 

2.476 


1.142 

1.146 


ros 

1.135 

ZBrNO 

114.4 

114.5 


ZBrON 

122.3 

Ac 

84326 

85500 


Ac 

96099 

Be 

3658 

3747 


Be 

2901 

Ce 

3506 

3586 


C c 

2816 

E 

0.130795 


BrBrO 

E 

0.105184 


0.017 



7? 

0.036 

fi 

1.00 




2.12 

f%xQ 

1.865 

1.843 


HlrBr 

2.510 

rsio 

1.865 

1.843 


HBrO 

1.690 

ZBrOBr 

112.9 

112.2 


ZBrBrO 

113.1 

A e 

32731 



A c 

17881 

Be 

1325 



Be 

1637 

C e 

1274 



C e 

1500 

E 

0.209868 

BrNC 

E 

0.148722 

V? 

0.015 




0.020 

u 

3.13 

2.94 


fi 

3.09 

r" 

rBrC 

1.810 

1.789 


HJrN 

1.800 


1.165 

1.158 


me 

1.183 

ZBrCN 

180 

180 


ZBrNC 

180 

A e 

4045 

4120 


Ac 

4332 

E 

0.243764 

CIBrO 

E 

0.221908 


0.016 



7? 

0.034 

u 

0.94 



fi 

2.26 

r L 

rcio 

1.730 



f ClBr 

2.322 


1.868 



TBiO 

1.682 

/ClOBr 111.7 



ZCIBrO 112.1 

Ac 

36563 



Ac 

18491 

B c 

2281 



Be 

3106 

Ce 

2147 



C t 

2659 

) E 

0.210970 




9? 

0.035 





M 

1.59 





CBrCl 

2.429 





rcio 

1.557 





ZBrCIO 115.1 





A c 

29245 





Be 

2192 





Ce 

2039 






« The energy is reported as — (E + 2648) for HOBr and HBrO, (E 
+ 2747) for FOBr and FBrO, -(£ + 2702) for BrNO and BrON, ~{E 
+ 5220) for BriO and BrBrO, -(£ + 2665) for BiCN and BrNCand 
_(£• + 3107 ) for ClOBr, CIBrO, and BrCIO. Obtained with the TZ2P 
hack set Bond lengths in angstroms and angles in degrees. Rotational 
constants are for the 79 Br isotope. » HOBr. substitution structure, dipole 
moment and vibrationally averaged rotational constants from ref 18. 
BrNO- vibrationally averaged structure from ref 10, dipole moment 
from refs 11 and 12, and vibrationally averaged rotational constants 
from ref 11. BrCN: vibrationally averaged structure and dipole 
moment from ref 19, vibrationally averaged rotational constant from 
ref 20. Br^O: vibrationally averaged structure from ref 45. c See rets 
17 and 41 for a detailed discussion of the 57T diagnostic. 


chlorine species, 1 ~ 5 it is expected that errors in the bond 
distances will be somewhat larger for hypervalent compounds 
than found for the molecules possessing normal Lewis dot 
structures. This is due to the fact that hypervalent species 
require larger one-particle basis sets to approach completeness 
relative to normal valent molecules. 

Agreement between CCSD(T)iTZ2P and experiment for the 
dipole moments of HOBr, BrNO, and BrCN is modest with 
the ab initio values being too high by between 0.11 0.26 D. 
The largest discrepancy exists for BrNO, but a close examination 
of the experimental determination shows that this value was 
guessed in a simulated spectrum. It thus seems likely that the 
true BrNO fi value is somewhat larger than 1.80 D, and further 
experimental work is needed. It is likely that the CCSD(T)/ 
TZ2P dipole moments are somewhat too high, as evidenced by 
the comparisons with experiment for HOBr and BrCN. This 
is due to one-particle basis set deficiencies and not to an 
inadequate treatment of electron correlation. Note that the 
experimental dipole moment for BrCN (3.13 D) is in very good 
agreement with the value obtained by Kello and Sadlej (3.19 
D) from high-level coupled-cluster calculations. Again, where 
no experimental data exist, the CCSD(T)fTZ2P dipole moments 
should aid in the analysis and interpretations of future experi- 
ments. 

It is interesting that all of the bromine triatomics studied here 
have rather sizable dipole moments— the largest is 4.00 D for 
HBrO and the smallest is 0.94 D for ClOBr. The dipole moment 
of the hypervalent species is always significantly larger than 
that determined for the analogous normal valent isomer, which 
is no doubt a reflection of the increased ionic character of the 
bonding (this is discussed in more detail later). 

B. Vibrational Frequencies and Quadratic Force Con- 
stants. CCSD(T)/TZ2P harmonic vibrational frequencies and 
IR intensities are presented in Table 2. Experimental data are 
given for comparison. Except for BrNO, the experimental data 
refer to fundamental vibrational frequencies and in some cases 
these are from matrix isolation experiments rather than gas- 
phase values. All vibrational data refer to the 79 Br isotope. 
Shifts relative to the 8, Br isotope are small, usually less than 1 
cm -1 . Of the vibrational modes included in Table 2, only the 
H-O stretch (str) and the H-Br str will have anharmonicities 
likely to be larger than 50 cm -1 . That is, the H-0 str 
anharmonicity is probably near 180 cm" 1 while the value for 
the H-Br str is less certain, but likely to be larger than 50 cm 
Taking into account the anharmonicity of the H—O str, then, 
the agreement between theory and experiment is very good for 
all modes which have been observed experimentally. This is 
true even for the BrBiO, CIBrO, and BrCIO hypervalent 
molecules. This good agreement for vibrational frequencies 
bolsters confidence in the reliability of the CCSD(T)/TZ2P 
predicted band centers where no experimental observations have 
been reported, and also in the reliability of the CCSD(T)/TZ2P 
equilibrium structures. 

In the study 27 where the Cl—O str of BrCIO and the Br— O 
str of CIBrO were observed, no bands from the ClOBr isomer 
were evident. This is in spite of the fact that ClOBr is the lowest 
energy isomer (see later discussion). The CCSD(TyTZ2P IR 
intensities given in Table 2 indicate that all of the fundamentals 
of ClOBr have rather small absorbances, and this explains why 
the ClOBr isomer was not observed. With regard to the IR 
intensities, another noteworthy result is that the N— O str in both 
BrNO and BrON has a very large IR intensity. It is hoped that 
the ab initio predictions in Table 2 will aid in future IR studies 
of bromine molecules. 
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TABLE 2: Harmonic Frequencies (cm 1 ) and IR In tensities (km/mol)* 

molecule/mode CCSD(T) expt* molecule/mode 


HOBr 

(O) (a) H— 0 str 

3807(73) 

i 02(0 ') bend 

1197(39) 

o^(a') 0— Br str 

608(7) 

FOBr 

<o\(a') F-0 str 

794 (33) 

(02(a) Br— 0 str 

605 (4) 

wi(a') bend 

297(1) 

BrNO 

<o\{a) N— 0 str 

1798(741) 

o ) 2 {a) bend 

545(57) 

0 )y(a) Br-N str 

268(34) 

Br 2 0 

a>j(a i) sym Br-0 str 

513(2) 

u>2{a\) bend 

180(0.2) 

a^to) asym Br-0 str 

613(0.1) 

BrCN 

o)i(o) C— N str 

2218(6) 

0)2 (0) Br— C str 

577(1) 

a>.i(;r) bend 

349(4) 

ClOBr 

o)\ (a) Cl-0 str 

658(1) 

0)2(0') Br— 0 str 

552(1) 

c 02(a) bend 

233(0.2) 

BrClO 

(0\ (a') Cl-0 str 

893(25) 

o)i{a'Y 

262(15) 

ujiia'y 

164(3) 


3610 

HBrO 

co](a) H— Br str 

1163 

co 2 (a') bend 

620 

o)3(a) Br-0 str 

1832 

FBrO 

0)] (a) Br-0 str 
c 02(a) F— Br str 
coy(a) bend 
BrON 

o)](a') N-0 str 

548 

o) 2(a) bend 

270 

coi(a') Br— O str 

526 

BrBrO 

oo](a') Br— O str 

180 

002(0') Br— Br str 

623 

g >3 (a) bend 

2198 

BrNC 

co](o) C— N str 

575 

C02 (a) Br-N str 

342 

o)i{n) bend 


CIBrO 

<o\(a ) Br-0 str 
Q)i{a) Cl— Br str 
cos(a') bend 

941 


CCSD(T) expt^ 


2292 (40) 

818(1) 

665(22) 


855(25) 

515(85) 

223(9) 


1827(1057) 

348(2) 

165(7) 


793(20) 

804 

215(16) 

153(3) 

236 

2089(124) 

561(0.2) 

194(0.2) 

2067 

808(25) 

301(34) 

182(4) 

819 


_ absotpuon intensmes in parentheses. All vibrational frequencies refer to the 79 Br isotope. 4 HOBr: fundamental frequencies from ref 21. 
BrNO. harmonic frequencies from ref 8. BrjO: stretching fundamental frequencies from ref 21 (matrix); bending fundamental from ref 45 BrBrO- 
fundtmenwl ^uencies from ref 2! (malnx). BrCN: fundamental frequencies from ref 20. BrNC: fundamental frequency from ref 21 (matrix). 
CIBrt) and BrClO: fundamenwlfrequencies from ref 27. c A potential energy distribution analysis shows that there is substantial mixine between 
the Br— Cl stretching and the BrCiO bending motions. 


Given the excellent agreement between experiment and theory 
for the vibrational frequencies, the CCSD(T)/TZ2P quadratic 
force fields are probably the most accurate available, and 
therefore these are presented in Table 3. The force constants 
should be useful in constructing force fields for larger bromine 
molecules. Units and internal coordinate definitions are given 
in a footnote to Table 3. Empirically derived force fields for 
BrNO and BrCN are given for comparison. The agreement 
between the CCSD(T)/TZ2P and empirical quadratic force fields 
is generally good. 

After submission of this manuscript, we became aware of a 
very recent experimental determination of the structure and 
quadratic force constants of B^O and the quadratic force 
constants of HOBr (see ref 45). For comparison, these new 
data have been included in Tables 1 -3. The agreement between 
the CCSD(T)/TZ2P results and the experimental determinations 
is very good. 

C. Bonding, Isomerization Energies, and Heats of For- 
mation. Mulliken populations of the various bromine triatomics 
are presented in Table 4. These are based on the CCSD(T) 
effective one-particle density matrics. As expected, the hyper- 
valent species exhibit a larger degree of ionic bonding as 
evidenced by larger atomic charges relative to the isomers with 
normal neutral Lewis dot structures. The ClOBr, CIBrO, and 
BrClO isomers provide a nice example of how hypervalent 
bonding involving Br is more pronounced than that involving 
Cl. That is, for CIBrO the Br is more positive than the Cl in 
BrClO, and the Cl partial charge is more negative in the former 
relative to the Br partial charge in the latter. Consistent with 
this, the Br partial charge in ClOBr is more positive than the 
Cl partial charge. It was pointed out in an earlier study on Cl 
hypervalent bonding 3 that the more electronegative the atom 
or moiety bonded to the hypervalent Cl atom, the more stable 
the hypervalent compound was relative to the normal valent 


isomer. Based on this conclusion and the above observations 
comparing Cl and Br partial charges, it is expected that Br 
hypervalent compounds will be more stable relative to the 
normal valent isomer as compared to their Cl analogs. Indeed 
this is found to be the case. 


Energies for the following isomerization reactions are given 
in Table 5. 

HOBr — HBrO + A£, 

(1) 

FOBr — FBrO + A E 2 

(2) 

ClOBr — CIBrO + A £ 3 

(3) 

Br,0 — * BrBrO + A£ 4 

(4) 

ClOBr — * BrClO + AE S 

(5) 

BrNO — BrON + A £ 6 

(6) 

BrCN — BrNC + A E 7 

(7) 


Only for reaction 3 are the MP2 energy differences consistently 
close to the CCSD( T) values (for the AN02-AN04 basis sets), 
but this is probably fortuitous since the CCSD energy differences 
are a few kcal/mol higher. The contribution of connected triple 
excitations is found to be small but significant for all of the 
reactions 1 —7. This suggests that the contribution from higher 
connected excitations should be very small. Examining the 
CCSD(T) energy differences, it is evident that the isomerization 
energies have nearly converged with respect to one -particle 
basis set improvements for all of the reactions. Based on the 
demonstrated convergence with respect to one- and n-particle 
basis set effects, it is expected that the CCSD(T)/AN04 energy 
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TABLE 3: Quadratic Force Constants Obtained at the 

CCSD{T)/TZ2P Level of Theory* 


molecule/ molecule/ 


constant 

CCSD(T) 

expC 

constant 

CCSD(T) 

HOBr 



HBrO 


fu 

8.09259 

8.085 

fu 

3.07313 

f\2 

-0.15490 

-0.168 

hi 

0.31116 

f\ 3 

0.04760 

0.103 

fn 

-0.15029 

/: 2 
fn 

3.14651 

3.363 

fn 

3.60873 

0.39706 

0.361 

fn 

0.20626 

fll 

FOBr 

0.74914 

0.738 

fn 

FBrO 

0.78771 

h i 

3.36303 


fn 

2.36098 

f] 2 

0.55998 


hi 

-0.05335 

fi 

0.26680 


fn 

-0.10933 

hi 

2.68299 


fn 

5.68835 

fiy 

0.19799 


fn 

-0.09264 

hi 

BrNO 

1.16227 


fn 

BrON 

0.63776 

fu 

14.60562 

15.25 ±0.04 

/n 

14.53247 

hi 

0.98918 

1.47 ±0.36 

hi 

-0.19092 

/.3 

0.21151 

0.11 ±0.20 

fn 

-0.03764 

fl2 

1.07106 

1.13 ±0.05 

fn 

0.40928 

hi 

0.09742 

0.10 ±0.02 

fn 

-0.00182 

hi 

BnO 

1.09265 

1.13 ± 0.02 

fn 

BrBrO 

0.43839 

/n 

2.53940 

2.812 

fn 

0.93948 

hi 

0.31766 

0.366 

hi 

-0.12918 

hi 

0.26614 

0.311 

fn 

-0.02792 

hi 

2.53940 

2.812 

f n 

4.89345 

hi 

0.26614 

0.311 

fn 

-0.05736 

hi 

BrCN 

1.08278 

1.026 

fn 

BrNC 

0.53196 

fu 

4.11647 

4.390 

fu 

3.87546 

hi 

0.03187 

-0.221 

hi 

-0.17974 

fl2 

17.48602 

17.194 

hi 

15.54058 

fll 

ClOBr 

0.32311 

0.304 

fn 

ClBrO 

0.10855 

fu 

2.80619 


fn 

1.24853 

f\ 1 

0.37960 


hi 

-0.08730 

hi 

0.29498 


fn 

-0.04480 

hi 

2.48549 


hi 

5.08879 

hi 

0.24138 


fn 

-0.06486 

hi 

BrCIO 

fu 

hi 

hi 

hi 

hi 

hi 

1.14917 

0.71243 

-0.24094 

-0.01511 

5.05397 

-0.03322 

0.56611 


fn 

0.58580 


a Units are aJ/A 2 , aJ/(A*rad), and aJ/rad 2 . Definition of the internal 
coordinates (including order) is as follows; HOBr: rno- r Bro< ZHOBr; 
HBrO: /’H B r. ''Bro, ZHBrO; FOBr: tfo, r B ro ZFOBr; FBrO: rrer, r Br o 
ZFBrO; BrNO: r N o, r B rN , ZBrNO; BrON: r N o, r Bl0 , ZBrON; Br 2 0: 
^BrO> ^BrOi ZBrOBr; BrBrO: r BrB r, r Br o ZBrBrO; BrCN: r Br c, fcN, ( 1 80 
- ZBrCN); BrNC: r BrN . rcN, (180 - ZBrNC); ClOBr: r C io, r Br o, 
ZClOBr; ClBrO: r OB r> r Br o, ZCIBrO; BrCIO: r BrC ,, r C \o , ZBrClO. 
b BrNO: ref 8. BrCN: ref 20. HOBr and Br 2 0: ref 45. 

differences should be accurate to better than ±0.8 kcal/mol. 
Adding the difference in zero-point vibrational energies (ZPVE) 
to the CCSD(T)/AN04 value, we obtain 56.3, -1 1.1, 1 1.3, 14.4, 
18.0, 23.7, and 35.9 kcal/mol for reactions 1-7, respectively. 
This is the reaction enthalpy at 0 K; for later purposes the change 
in the reaction enthalpies due to thermal population of vibra- 
tional energy levels at 298.15 K is also presented in Table 5. It 
is thus apparent that the isomers having a normal neutral Lewis 
dot structure are the lowest energy species in every case except 
for FOBr/FBrO. This is similar to the situation for the Cl 
analogs (see refs 1-5 and references therein), but as indicated 
previously, bromine hypervalent isomers are more stable than 
their chlorine analogs. Again the ClOBr, ClBrO, and BrCIO 
isomers provide a useful example. While ClOBr is the lowest 


TABLE 4: Atomic Charge Distribution Based on Mulliken 

Population Analyses* 



HOBr 

HBrO ClOBr 

ClBrO 

BrCIO 

H/Cl 

±0.33 

±0.02 ±0.14 

-0.15 

±0.50 

O 

-0.52 

-0.62 -039 

-0.52 

-0.43 

Br 

±0.19 

±0.60 ±0.25 

±0.68 

-0.07 


FOBr 

FBrO 

Br 2 0 

BrBrO 1 ' 

F/Br 

-0.12 

-0.39 

±0.24 

-0.09 

O 

-0.16 

-0.55 

-0.48 

-0.51 

Br 

±0.29 

±0.94 

±0.24 

±0.60 


BrNO 

BrON 

BrCN 

BrNC 

O/C 

-0.14 

-0.06 

-0.28 

±0.12 

N 

±0.33 

±0.42 

-0.14 

-0.49 

Br 

-0.19 

-0.36 

±0.42 

±0.37 

" Based on the CCSD(T)/TZ2P effective one-particle density ma- 

trices. * 

The third row refers to the central Br atom. 


TABLE 5: Isomerization Energies (kcal/mol) 0 



A£, 

A £: A£ 3 ■ 

A E a A£ 5 

A£ 6 A Ei 


MP2/AN01 

64.3 

-11.9 

14.7 

19.7 

21.8 

19.9 

41.9 

MP2/AN02 

61.7 

-15.5 

12.6 

17.7 

20.4 

21.3 

41.4 

MP2/AN03 

60.4 

-17.3 

11.7 

17.0 

19.7 

21.5 

41.3 

MP2/AN04 

60.5 

-17.4 

11.5 

16.7 

19.5 

21.6 

41.3 

CCSD/ANOl 

62.7 

-4.0 

17.8 

20.8 

25.0 

27.3 

37.1 

CCSD/AN02 

60.7 

-7.0 

16.1 

19.5 

24.2 

28.7 

36.7 

CCSD/AN03 

59.7 

-8.5 

15.6 

19.1 

24.0 

29.0 

36.6 

CCSD/AN04 

60.2 

-8.5 

15.5 

19.0 

23.9 

29.2 

36.7 

CCSD(T)/TZ2P 

61.5 

-4.7 

13.7 

16.1 

20.6 

22.6 

38.4 

CCSD(T)/AN01 

61.6 

-6.3 

13.9 

16.6 

19.5 

21.9 

37.3 

CCSD(T)/AN02 

59.7 

-9.2 

12.3 

15.2 

18.7 

23.1 

36.9 

CCSD(T)/AN03 

58.4 

-10.9 

11.6 

14.8 

18.4 

23.6 

36.8 

CCSD(T)/AN04 

58.9 

-11.0 

11.5 

14.6 

18.2 

24.1 

36.8 

AZPVE'’ 

-2.6 

-0.1 

-0.2 

-0.2 

-0.2 

-0.4 

-0.7 

AVT 

0.0 

0.1 

0.2 

0.2 

0.2 

0.2 

0.3 


See eqs 1-7 for definition of the A £ quantities. Energies 
determined at CCSD(T)/TZ2P equilibrium geometries. Zero-point 
vibrational energies not included-see text for energy differences where 
these are included. b Zero-point vibrational energies determined using 
CCSD(T)/TZ2P harmonic frequencies. c Energy changes due to thermal 
population of vibrational energy levels at 298.15 K. 

energy isomer, comparison of reaction energies A £3 and A£s 
(above) shows that ClBrO is about 6.7 kcal/mol more stable 
than BrCIO. Further, comparison of the isomerization energies 
of the analogous Cl compounds supports this assertion (e.g., 
FCIO is predicted 3 to be more stable than FOCI by only 6.0 
kcal/mol rather than the 11.1 kcal/mol found for FBrO/FOBr). 

Comparison of A E ? through A£V shows that for XBrO 
compounds, the more electronegative the X species the more 
stable the hypervalent compounds, similar to chlorine hyper- 
valent compounds. The accurate isomerization energies deter- 
mined from the data in Table 5 should be useful in assessing 
the stratospheric significance of bromine compounds as well 
as determining accurate heats of formation. However, in order 
to evaluate accurate heats of formation from this data it is first 
necessary to place some of these bromine compounds on an 
absolute scale. That is, it is necessary to determine accurate 
heats of formation for a few bromine compounds first. In 
previous studies on fluorine and chlorine compounds, our 
approach has been to compute accurate isodesmic and/or 
homodesmic reaction energies and then use these together with 
accurate experimental heats of formation of some compounds 
in order to determine heats of formation of new compounds. 
This approach has the advantage that the energy differences 
determined from the ab initio calculations involve a high degree 
of cancellation of errors, and therefore it is relatively easy to 
approach one- and n-particle space convergence (note that the 
cancellation of errors is certainly better for isodesmic compared 
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to homodesmic reactions, but with CCSD(T) and large ANO 
basis sets residual errors are generally small even for homodes- 
mic reactions). In the case of chlorine and especially bromine, 
there is also the big advantage that since only closed-shell 
molecules are used in the isodesmic and homodesmic reactions, 
it is not necessary to include explicitly spin— orbit effects. The 
difficulty with this approach is that it requires knowledge of an 
accurate experimental heat of formation for at least one closed- 
shell bromine. As indicated previously, there is some debate 
currently 22 " 24 concerning the heat of formation of HOBr, with 
estimates ranging from —9 to —19 kcal/mol (298 K), although 
the more recent determinations range from —9 to —14 kcal/ 
mol. This range is still quite large, but probably the A//°f ?98 
(Atf°f,o) value determined by Ruscic and Berkowitz, 23 —13.4 
± 0.4 kcal/mol (—10.9 ± 0.4 kcal/mol; these are actually lower 
bounds), is the most reliable since the auxiliary thermochemical 
data used in this study is the most reliable (relative to that used 
in refs 22 and 24). According to ref 14, the heat of formation 
of BrNO is 19.6 ± 0.2 kcal/mol (298 K), which has been 
obtained from experiments carried out in the 1920s and 1930s. 
This is a very small uncertainty, although there are many 
examples in which the heats of formation of similar fluorine 
and chlorine compounds were thought to possess very small 
uncertainties and this was incorrect (e.g., see refs 42 and 43). 
We may examine the consistency of the HOBr and BrNO 
experimental heats of formation and also determine heats of 
formation for the other bromine compounds included in this 
study by computing the reaction energies of the following 
homodesmic (reactions 8 and 12) and isodesmic (reactions 
9—11) reactions, and these are presented in Table 6. 

BrNO + H 2 0 — HNO + HOBr + A £ 8 (8) 

HOBr + HOBr — H 2 0 + Br 2 0 4- A£ 9 (9) 

Br 2 0 + HOF — HOBr + FOBr + A£ 10 (10) 

HOBr + HOC1 — ClOBr + H 2 0 + A E u (11) 

BrCN + HOCl — C1CN + HOBr + A£ 12 (12) 

A quick examination of the results in Table 6 shows that the 
reaction energies A E% through A En are very nearly converged 
with respect to one- and n-particle basis set improvements. 
Therefore, taking the CCSD(T)/AN04 A E values as our best 
computed result and correcting these with the contribution from 
zero-point energies, we obtain 46.5, —2.3, —2.6, —2.4, and —8.3 
kcal/mol as the best estimates for AE% through AE 12 , respectively 
(0 K). Based on the exhibited convergence and the established 
reliability of the CCSD(T) method, these values are expected 
to be accurate to better than ±1.0 kcal/mol. Combining AEs 
with experimental heats of formation for BrNO and H 2 O , 14 and 
a new determination 44 for the heat of formation of HNO (26.7 
and 26.0 kcal/mol for A//°f,o and AW°f.298» respectively), the 
heat of formation of HOBr is computed to be —15.5 kcal/mol 
(0 K) or —18.0 kcal/mol (298 K). The Atf°f.298 value is not in 
very good agreement with the determination of Ruscic and 
Berkowitz, being 4.6 kcal/mol lower. This is especially bad 
agreement considering that the Ruscic and Berkowitz value is 
a lower bound. Hence it is reasonable to conclude that the 
experimental A//°f ,298 value for BrNO is in error and instead 
we use the experimental heats of formation of HOBr 23 and 
H 2 O , 14 the recent determination for HNO , 44 and the best estimate 
for A£g to arrive at a new A//°f ,298 value for BrNO, 26.4 kcal/ 
mol. It is difficult to assign an uncertainty to the present 


TABLE 6: Homodesmic and isodesmic reaction energies 

(kcal/mol) 0 



A £9 

A£io 

A£ n 

A£ 12 

MP2/ANO 1 

53.3 -3.5 

-0.9 

-3.2 

-8.5 

MP2/AN02 

50.8 -3.4 

-0.9 

-3.1 

-9.0 

MP2/AN03 

51.2 -3.4 

-1.2 

-3.2 

-9.0 

MP2/AN04 

51.1 -3.4 

-1.1 

-3.2 

-9.2 

CCSD/ANOl 

45.5 -0.9 

-2.0 

-1.0 

-7.6 

CCSD/AN02 

42.8 -0.7 

-2.0 

-0.8 

-8.2 

CCSD/AN03 

43.2 -0.7 

-2.4 

-0.7 

-8.1 

CCSD/AN04 

43.2 -0.6 

-2.3 

-0.6 

-8.4 

CCSD(T)/TZ2P 

49.3 -1.2 

-2.8 

-1.3 

-6.4 

CCSD(T)/AN01 

49.2 -1.8 

-2.3 

-1.9 

-7.5 

CCSD(T)/AN02 

46.8 -1.7 

-2.2 

-1.8 

-8.1 

CCSD(T)/AN03 

47.1 -1.7 

-2.6 

-1.7 

-8.0 

CCSD(T)/AN04 

47.1 -1.6 

-2.5 

-1.7 

-8.3 

AZPVE* 

-0.6 -0.7 

-0.1 

-0.7 

0.1 

AVT C 

-0.3 0.4 

-0.1 

0.4 

0.0 

See eqs 8—12 for definition of the A E quantities. Zero-point 
vibrational energies not included— see text for energy differences where 

these are included. 

b Zero-point vibrational energies determined using 

CCSD(T)/TZ2P harmonic frequencies. c 

Energy changes due to thermal 

population of vibrational energy levels at 298.15 K. 



TABLE 7: Heats of Formation (kcal/mol) 0 




Atf° f .o 


A//°f.29S 15 


HBrO 

+45.4 


+42.9 


FOBr 

+21.4 


+ 19.2 


FBrO 

+ 10.3 


+8.2 


ClOBr 

+26.7 


+24.6 


CIBrO 

+38.0 


+36.1 


BrClO 

+44.7 


+42.8 


Br 2 0 

+33.0 


+29.1 


BrBrO 

+47.4 


+43.7 


BrNO 

+26.4 


+24.2 


BrON 

+50.1 


+48.1 


BrCN 

+47.2 


+45.6 


BrNC 

+83.3 


+82.0 



0 Based on the data in Tables 5 and 6 together with some 
experimental and theoretical heats of formation— see text. 


computed heat of formation— the estimated uncertainties of all 
of the individual contributing components are rather small, and 
since these are unrelated the overall uncertainty should be less 
than their sum. We believe that the uncertainties assigned to 
the computed reaction energies are conservative, so if the 
A//° f ,298 value for BrNO determined in this work is in error by 
several kcal/mol, then the largest error will probably arise from 
the experimental heat of formation of HOBr. 

In any case, assuming that the experimental heat of formation 
for HOBr 23 and the computed heat of formation for BrNO are 
accurate, the data contained in Tables 5 and 6, together with 
experimental heats of formation for H 2 O, HOF, 43 HOCl, 14 and 
C1CN, 14 may be used to determine heats of formation for the 
remaining bromine compounds included in this work. These 
are collected in Table 7. It is noteworthy that Orlando and 
Burkholder 24 were able to deduce that AJ/KB^O) — 2A Hr 
(HOBr) = 55.5 kcal/mol (assuming that AS ^ 0) based on an 
equilibrium between HOBr, Br?0, and H 2 O being reached. In 
the present work we obtain 55.9 kcal/mol for this quantity (from 
reaction 9 and the experimental AJfV^s^O) value), which 
strongly supports the accuracy of the present ab initio calcula- 
tions. Should an unquestionably accurate experimental heat of 
formation of any of the bromine triatomics included in this study 
be obtained, the heats of formation of all of the other bromine 
triatomics may be derived from the data in Tables 5 and 6 
together with the previously noted experimental heats of 
formation of H 2 O, HOF, HOCl, and C1CN and the new 
theoretical heat of formation of HNO. 44 
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Conclusions 

The CCSD(T)/HTZ2P level of theory has been used to 
determine accurate equilibrium geometries, harmonic frequen- 
cies, dipole moments, and quadratic force fields for HOBr, 
HBrO, FOBr, FBiO, ClOBr, CIBrO, BrCIO, Br 2 0, BrBrO, 
BrNO, BrON, BrCN, and BrNC. The CCSD(T)/TZ2P geom- 
etries and vibrational frequencies are in very good agreement 
with the available experimental data, although the CCSD(T) 
dipole moment values are consistently smaller than experiment. 
This is most likely due to errors in both the theoretical treatment 
(i.e., one-particle basis set deficiencies) and the experimental 
analyses. An explanation is provided for the fact that ClOBr 
was not observed in matrix isolation experiments in which 
CIBrO and BrCIO were both detected. That is, the IR intensities 
of the stretching modes in ClOBr are very small. It is expected 
that the ab initio structures, vibrational frequencies, and 
quadratic force fields should be useful in the analysis of future 
experimental studies. 

Accurate isomerization energies are computed at the CCSD- 
(T) level of theory using large ANO basis sets that include 
g-type functions. With one exception, the isomer possessing a 
normal neutral Lewis dot structure is the lowest energy species. 
In the one exception FBrO is found to be 11.1 kcal/mol lower 
in energy than FOBr (0 K). Using a set of homodesmic and 
isodesmic reaction energies, heats of formation of all of the Br 
triatomics have been computed. The accuracy of the ab initio 
calculations is demonstrated by the excellent agreement with a 
recent experimental study 24 for the determination of the A/f f - 
(Br 2 0) - 2A//f{HOBr) quantity. The accuracy of the proposed 
heats of formation, however, ultimately depends on the reliability 
of the experimental heat of formation of HOBr. 
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